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Abstract 

The collimation C of a hadronic event in the e + e~-annihilation is de- 
fined as the average of cos#, C =< cos# >, where is the angle of each 
hadron measured from the thrust axis, and the average is over all the 
hadrons produced in an event. It is an infrared-stable event-shape param- 
eter. 1 — C, the difference between the unity and the average collimation at 
a given energy, is proportional to the anomalous dimension of the hadron 
multiplicity at the leading order in MLLA. Its next-to-leading order cor- 
rections are calculated. 
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When an energetic coloured particle (parton) is produced at a high-energy 
collision, many hadrons are produced near the direction of the energetic parton 
(jet phenomenon). The shape of a jet provides us significant information on the 
dynamics underlying the multiple hadroproduction. 

In order to make quantitative predictions by the perturbative QCD, and to 
test them against the experiments, the quantity which characterises the event 
shape (event-shape parameter) has to be infrared(IR)-stable at asymptotically 
high energies: i.e. insensitive to the IR cutoff under which the QCD interaction 
becomes non-perturbative. Otherwise, the quantity would depend on the non- 
perturbative hadronisation process, on which our theoretical knowledge is very 
limited. 

The thrust T in the e + e~-annihilation is an example of the IR-stable event- 
shape parameter. The thrust distribution is known to be IR-stable even near 
T = 1 (collinear limit) §0. 

In this article, we introduce another IR-stable event-shape parameter, the 
collimation C of a hadronic event in the e + e~-annihilation. It is defined event by 
event by 

C=<cos9>=— (1) 

n 

where 9{ is the angle of each hadron measured from the thrust axis. The average 
is over all the n hadrons in an event. C = 1 in the limit that all the hadrons are 
in parallel, while C — 1/ 2 in the limit that the angular distribution of hadrons 
is spherically symmetric. 

Let us analyse 1 — C, the difference between the unity and the average colli- 
mation at a given energy. It is defined by 

2/^2(1-77)5 



C = l-cosfl= ^ _ ^ (2) 
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where n is the average hadron multiplicity, and dn/ dr] its angular density. d<7( 1P ;n / dr/ 
is the one-particle-inclusive cross section of a hadron, and <7o the total hadronic 
cross section. Let us call 1 — C the "jet-width"Q. 

At asymptotically high energies, the average multiplicity and its angular den- 
sity are determined by the successive soft-gluon emission, and can be evaluated 
by the resummation of the associated double-logarithms to all orders in a s (the 
modified leading-log approximation. MLLAjlj || || ||). 

As we shall see below, the jet-width is proportional to the anomalous dimen- 
sion 7 of the multiplicity at the leading order in MLLA: 

1 - C = (2 In 2)7 + (higher order) (3) 



2 The jet-width is determined by the average angular distribution of hadrons at a given 
energy. C, therefore, is not exactly identical to the average of the value of the collimation C 
measured for each hadronic event. 
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where 7 is defined by the dependence of n on the centre-of-mass energy W of the 

event _ 

dlnn . . 

7 = din ' ^ ^ 

We take the QCD mass-scale Aqcd as the unit of energy for the notational 

simplicity. 

If the local parton-hadron duality (LPHDH p]]) holds for the particle flow, the 
modification of the jet-width by the hadronisation process disappears at high en- 
ergies. The measurement of the jet-width thus provides us a test of the theoretical 
framework of MLLA+LPHD. 

The anomalous dimension 7, which is proportional to v/aj, is the basic quan- 
tity in MLLA. It is the expansion parameter in the resummed perturbation the- 
ory. 7 can be measured directly by comparing the multiplicity data at different 
energies (with the use of its definition (4)). The result, however, remains some- 
what ambiguous because of the uncertainty in the systematic errors in different 
experiments. One of the advantages of using (3) is that the value of 7 is deter- 
mined at each experiment. (The normalisation constant of n, which cannot be 
determined by the perturbation theory, is cancelled between the numerator and 
the denominator on the rhs of (2).) We thus have a less ambiguous comparison 
of the prediction with the experimental data. 

The determination of 7 by measuring the jet- width 1 — C, however, is indirect: 
The proportionality is only at the leading order. It is, therefore, important to 
evaluate the higher order corrections. 

Roughly speaking, the jet-width is the ratio of the multiplicity in the large- 
angle region to the total multiplicity. The multiplicity at large angles is also 
interesting because the anglular ordering fj || [|T0l , on which the shower Monte 
Carlo simulation programs EI pi are based, needs modification at large angles 



(so-called the "dipole corrections" P| |TT|). 

The average multiplicity in the e + e~-annihilation is written in terms of the 
multiplicity in a gluon-subjet, M g (k\ 
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where Ca = 3 is the gluon colour-charge, Cp = 4/3 the quark colour- charge, 
and Nf the number of active quark flavours, k is the momentum of the gluon 
emitted from the quark- ant iquark pair (e + e~ — > qq + g), and Q\ its polar angle 



3 



(the direction of the antiquark is chosen as the — ^-direction). da/(dr]idz) is the 
one-gluon emission cross section at 0(a s ): 

1 da C F a s 2(1 -z)(l- z^) 2 + (1 - zfrjjz 2 + z 2 {l - Vl ) 2 



0(a s ) 



7T ZT]i(l - ?7i)(l - 2zT]i + Z 2 T]i) 



The normalisation constant A on the rhs of (6) is not determined by the pertur- 
bation theory, and is cancelled between the numerator and the denominator in 
the expression for the jet-width (2). For simplicity, we put A = 1 below. 

The integration over the gluon momentum on the rhs of (5) (with (6) and (7)) 
gives at the next-to-leading order 



n 



^f 27l 1/2+Bo7g | 
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B = --ln2. 
2 

The jet-width (2) is given at the leading order in MLLA by 

o ri r 1 Hit 

1-C=-/ d m 2(l- m ) dz--—M 9 (kl), (9) 
n Ji/2 7 Go dr^id^ 

In (9), we identified the direction of the registered hadron to the direction of 
the gluon-subjet[|J. It was shown in JTTJ that the correction due to the difference 



between the two directions vanishes at the next-to-leading order. 

In (9), as well as in (5), the evaluation is made for the forward emission only 
(rj > 1/2), and the result is doubled. Let us assume that the gluon is less energetic 
than the antiquark. (The case that the gluon is the most energetic among the 
initial three partons shall be analysed later.) Then the condition rj > 1/2 implies 
that the antiquark is the most energetic and that the thrust axis is in the direction 
of the antiquark. The upper bound for the gluon energy- fraction z = 2k/ W is 
now 

2771-1 1 3. 

z < (for-<r 7l < i ), 

z < 1 - V 1 - Vl (for - < 77! < 1 ) . (10) 
Vi 4 

Evaluation of (9) at the next-to-leading order (with (6), (7) and (10)) gives 

(9 ) = ^L 2 fi- J- 1 ' 2 U , B ^ }^gVH, ( u; 

C A n [ y 7ov /ln(^ 2 /12)y ^J\n(W 2 / 12) J ^J\n(W 2 / 12) 

Bl = _H_9 lll2 + | lll3 4 + isp(i) + is P (5)-sp(I)-Sp(I 

Sp(x) = / —In- . 

Jo y 1-y 
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With the use of (8) for n, we obtain at the next-to-leading order 
1 - cos 9[ = 2 In 2 {1 + (- In 3 + £0)7} 7 + 25i7 2 , 
where 7, defined by (4), is 



7 = 



7o 



+ 



7i 



ln(W 2 /12) ln(l^ 2 /12) 



(12) 



(13) 



So far, we have considered the one-gluon emission (from the initial quark- 
antiquark pair), which produces a gluon-subjet to which the registered hadron 
belongs. When another hard gluon is emitted, its recoil may change the direction 
of the parent quark. For such additional hard-gluon emission, its amplitude has 
an extra factor of a s . Because the jet-width is of O^^/al), the 0(a s ) contribution 
is at the next-to-leading order. The angle of the gluon-subjet is expressed as the 
vector-sum of the recoil angle 9 2 of the parent quark and the gluon emission angle 
61 from the quark: 9 = 61 + 9 2 . The jet-width is thus 



1 - C = 1 - cos9 = 1 - COS01 + 1 - cos# 2 + (1 - COS0i)(l - cosfl 2 ) • (14) 

The third term on the rhs of (14) is 0(7 3 ). At the next-to- leading order, therefore, 
we only need to evaluate the first and the second terms separately. In terms of 
the recoil angle 9 2 and the energy fraction z\ = 2P 1 / W (Pi. the energy of the 
quark), the one-gluon emission cross section at 0(a s ) reads 



da 



d^d^i 



o(a 3 ) 



V2 



x 



2n 1 - 772 
1 + cos 9 2 



1 - 2zxr] 2 + z\r} 2 
+ zi(l - zi)(3?72 - 



+ l-z 1 



(15) 



By the use of (15), the average value of 1 — cos 9 2 is calculated to be 



cos 



9 2 = — I' d^ 2 2(l - 772) / cbi 

<7n Jl/2 J 
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Bo = 




di] 2 dz 1 
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3 4 
7 In - + In 2 In - 

2 3 



5s P (i) + isp(i) + Sp^ 



where the integration over the quark momentum is done under the condition that 
the antiquark is the most energetic among the three initial partons. 

Finally, let us analyse the case that the gluon is the most energetic among the 
initial three partons. The thrust axis is now identical to the direction of the hard 
gluon. The soft-gluon emission parallel to the gluon would contribute to n at 
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0(a s ), but does not contribute to the jet-width. The soft-gluon emission parallel 
to the (anti) quark, on the other hand, gives a next-to- leading contribution to the 
jet-width. At this order, therefore, the direction of the registered hadron can be 
identified to the direction of the (anti) quark. In terms of the polar angle 9 3 of 
the quark (the gluon is now in the — 2-direction), the one-gluon emission cross 
section at 0(a s ) reads 



da 



dr] 3 dz 1 



0(a s ) 



C F a s ( 2:1(1 - zi) 1- zi 

2vr 1 - 224773 + zfr] 3 \ 1 - 21773 r] 3 

1 + COS #3 

Vs = g • 

By the use of (17), the contribution to the jet- width is calculated to be 



If 1 r da 

1- cos 6> 3 = — / d77 3 2(l - 773) / dzi 
an J1/2 J 



a J1/2 J dr] 3 dz 1 



o(a s ) 71 



C F a s 
-D3 



B 3 = (-— -61n2 + ^ln3 + 211n 2 2 + 21n 2 3- 161n21n3 
12 4 

-4Sp(l) + 6Sp Q - 2Sp (1)](18) 

From (12), (16) and (18), we obtain 



1 — cos 9 = 1 — cos ^1 + 1 — cos 9 2 + 1 — cos 6* 3 

= (2 In 2)7 + {-2 In 2 In 3 + 2(5 In 2 + 5x) } 7 2 + (£ 2 + 509) 



2 , ^Vo.s 



7T 



Numerically, the factors in (19) are 2(£ ln2 + -B 1 ) ~ -0.086, B 2 w 0.432 and 
i? 3 0.128. The origin of the largest part of the next-to-leading order correction, 
— 2(hi21n3)7 2 (in the second term on the rhs of (19)), is the upper bound of the 
gluon energy in evaluating (9). In fact, as the gluon energy fraction z increases 
beyond the limit of (10), the event-shape becomes more collimated (the thrust 
axis is now in the direction of the gluon). 

The prediction for the jet-width is listed in Table 1. In the prediction, (9) 
and (7) are evaluated by performing the numerical integration, rather than by 
the use of the expansion in 7 (12). The corresponding value of the anomalous 
dimension is given by (13). The value of (one-loop) QCD mass scale is chosen 
to be Aqcd = 0.15GeV. The dependence on the cm energy W (hence on Aqcd) 
is rather weak over wide range of W. This is because it decreases proportionally 
to 1/Vln W 2 . Indeed, we cannot use the prediction for the determination of 
the universal (but scheme- dependent) QCD mass scale (such as A^r) until the 
next-to-next order correction is calculated. 
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Table 1 

The prediction for the jet-width at various energies, and the corresponding 
value of the anomalous dimension. 



E 


1- < C> 


7(^712) 


58GeV 


0.251 


0.236 


92GeV 


0.244 


0.228 


lOOOGeV 


0.216 


0.195 
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